Laboratory Tests of Low Density Astrophysical Equations of State 
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Clustering in low density nuclear matter has been investigated using the NIMROD multi-detector 
at Texas A&M University. Thermal coalescence modes were employed to extract densities, p, and 
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temperatures, T, for evolving systems formed in collisions of 47 A MeV Ar + 



Sn and 



b4 Zn + x ^Sn , ^ 4 Sn. The yields of d, t, *Re, and 4 He have been determined at p = 0.002 to 0.03 
nucleons/fm 3 and T= 5 to 11 MeV. The experimentally derived equilibrium constants for a particle 
production are compared with those predicted by a number of astrophysical equations of state. The 
data provide important new constraints on the model calculations. 
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I. INTRODUCTION 



Reliable understanding of the nuclear equation of state, 
EOS, over a wide range of densities and temperatures is 
crucial to both nuclear science and to our understand- 
ing of stellar evolution and supernovae [l[ . In the latter 
context it is well known that a valid treatment of the 
correlations and clusterization in low density matter is 
a vital ingredient of astrophysical models. To meet the 
need for the nuclear input, some extensive well known 
calculations and existing tabulations, based on varying 
effective interactions, were developed and have served as 
standard input for a wide variety of astrophysical simu- 
lations [2h4J • More recently some new approaches have 
led to a number of new predictions [a-llul]. While all of 
the models in use lead to strong alpha clustering of the 
matter at low nucleon densities, p, and temperatures, T, 
they differ significantly in their quantitative predictions, 
usually tabulated as alpha mass fractions, at specified 
T, p and proton mass fraction, Yp. See references @, IH 
for example. For a given Yp the differences in the abso- 
lute values reflect differences in the effective interactions 
inherent in the chosen nuclear equation of state, differ- 
ences in the formulation and approximations employed in 
the models, and differences in the treatment of possible 
competing species. We present here results of an exper- 



imental determination of clustering yields in low density 
nuclear matter and use these results to make a direct 
test of the different models, focusing on alpha clustering. 
Our test observable is not the alpha mass fraction but 
rather the equilibrium constant for alpha particle clus- 
tering. The model derived equilibrium constants should 
be nearly independent of proton fraction and the choice 
of competing species assumed in a particular model. 



EXPERIMENTAL TECHNIQUES 

We reported in Refs. [lH and that measurements 
of nucleon and light cluster emission from the participant 
matter which is produced in near Fermi energy heavy ion 
collisions could be employed to probe the EOS at low den- 
sities and moderate temperatures where clustering is im- 
portant. Our data demonstrated a large degree of alpha 
clustering for densities at and below ~ 0.05 times normal 
nuclear density and temperatures of 4 to 10 MeV. Using 
these data we derived experimental symmetry free ener- 
gies in low density nuclear matter ll|, ll2| ■ That analysis 



employed the isoscaling technique which compares yields 
for two systems with similar temperatures but different 
N/Z ratios to determine the differences in chemical po- 
tentials and symmetry energy 3, HI- 
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The NIMROD 4tt multi-detector at Texas A&M Uni- 
versity has now been used to extend our clustered mat- 
ter measurements to higher densities. Cluster produc- 
tion in collisions of 47A MeV 40 Ar with 112 ' 124 Sn and 
64 Zn with 112 < 124 Sn was studied. NIMROD consists of a 
1G6 segment charged particle array set inside a neutron 
ball |l5[. The charged particle array is arranged in 12 
rings of Si-Csl telescopes or single Csl detectors concen- 
tric around the beam axis. The Csl detectors are 1-10 cm 
thick Tl doped crystals read by photomultiplier tubes. A 
pulse shape discrimination method is employed to iden- 
tify light particles in the Csl detectors. For this exper- 
iment each of the forward rings included two segments 
having two Si detectors (150 and 500 /im thick) in front 
of the Csl detectors (super telescopes) and three having 
one Si detector (300 (im thick). Each super telescope was 
further divided into two sections. Neutron multiplicity 
was measured with the 4-7T neutron detector surrounding 
the charged particle array. This detector is a neutron 
calorimeter filled with Gd doped pseudocumene. Ther- 
malization and capture of emitted neutrons in the ball 
leads to scintillation which is observed with phototubes 
providing event by event determinations of neutron mul- 
tiplicity. Further details on the detection system, energy 
calibrations and neutron ball efficiency may be found in 
Ref. [16J. The combined neutron and charged particle 
multiplicities were employed to select the most violent 
events for subsequent analysis. 



ANALYSIS 

The dynamics of the collision process allow us to probe 
the nature of the intermediate velocity "nucleon-nucleon" 
17-19j]. Measurement of emission cross 



emission source 
sections of nucleons and light clusters together with suit- 
able application of a coalescence ansatz [20| provides the 
means to probe the properties and evolution of the in- 
teraction region. The techniques used have been de- 
tailed in several previous publications [0,12,17- 19| and 
are described briefly below. A notable difference from 
Refs. ll|, ll2[ is the method of density extraction. This is 
discussed more extensively in the following. We empha- 
size that the event election is on the more violent colli- 
sions. Cross section weighting favors mid-range impact 
parameters. 

An initial estimation of emission multiplicities at each 
stage of the reaction was made by fitting the observed 
light particle spectra assuming contributions from three 
sources, a projectile-like fragment (PLF) source, an inter- 
mediate velocity (IV) source, and a target-like fragment 
(TLF) source. A reasonable reproduction of the observed 
spectra is achieved. Except for the most forward detec- 
tor rings the data are dominated by particles associated 
with the IV and TLF sources. The IV source veloci- 
ties are very close to 50% of the beam velocity as seen in 



many other studies [17H19j . The observed spectral slopes 

reflect the evolution dynamics of the source [a nana. 

For further analysis, this IV source is most easily sam- 
pled at the intermediate angles where contributions from 
the other sources are minimized. For the analysis of the 
evolution of the source we have selected the data in ring 
9 of the NIMROD detector. This ring covered an angu- 
lar range in the laboratory of 38° to 52°. Inspection of 
invariant velocity plots constructed for each ejectile and 
each system, as well as of the results of the three-source 
fit analyses indicate that this selection of angular range 
minimizes contributions from secondary evaporative de- 
cay of projectile like or target like sources. We treat the 
IV source as a nascent fireball created in the participant 
interaction zone. 

The expansion and cooling of this zone leads to a cor- 
related evolution of density and temperature which we 
probe using particle and cluster observables, yield, en- 
ergy and angle. As in the previous work II lL 1 1 21 ] we have 
employed double isotope yield ratios 23|,|24( to character- 
ize the temperature at a particular emission time. Model 
studies comparing Albergo model (23j temperatures and 
densities to the known input values have shown the dou- 
ble isotope ratio temperatures to be relatively robust in 
this density range [251 ] . However the densities extracted 
using the Albergo model are useful only at the very lowest 
densities [25| . Consequently, in this study we have em- 
ployed a different means of density extraction, the ther- 
mal coalescence model of Mekjian I?], 20j ] . 

To determine the coalescence parameter Pq, the radius 
in momentum space, from our data we have followed the 
Coulomb corrected coalescence model formalism of Awes 



et aV 26] and previously employed by us in Ref. |17|. 
In the laboratory frame the derived relationship between 
the observed cluster and proton differential cross sections 
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where the double differential multiplicity for a cluster of 
mass number A containing Z protons and N neutrons 
and having a Coulomb-corrected energy Ea, is related 
to the proton double differential multiplicity at the same 
Coulomb corrected energy per nucleon, E — Ec, where 
Ec is the Coulomb barrier for proton emission. R np is 
the neutron to proton ratio. 

A strict quantitative application of the coalescence 
model requires knowledge of cluster, neutron and proton 
differential cross sections with proper absolute normal- 
izations. In this work absolute measured multiplicities 
for the selected violent events are employed. The neu- 
tron spectra are not measured. However, since within 
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FIG. 1: Temperatures and densities sampled by the expand- 
ing IV source. The shaded band indicates the uncertainity in 
density. See text. 



the framework of the coalescence model the yield ratios 
of two isotopes which differ by one neutron are deter- 
mined by their binding energies and the n/p ratio in the 
coalescence volume, we have used the observed triton to 
3 He yield ratio to derive the n/p ratio used in this anal- 
ysis. 

In the Mekjian model thermal and chemical equilib- 
rium determines coalescence yields of all species. Under 
these assumptions there is a direct relationship between 
the derived radius in momentum space and the volume 
of the emitting system. In terms of the Po derived from 
Eq. ([T]) and assuming a spherical source, 
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where h is Plancks constant and Z, N, and A are the 
same as in Eq. ([T]) . E is the binding energy, s the spin of 
the emitted cluster, and T is the temperature. Thus the 
radius, and therefore the volume, can be derived from the 
observed Po and temperature values assuming a spherical 
shape. We note that this volume is a free volume. 

Because our goal was to derive information on the den- 
sity and temperature evolution of the emitting system, 
our analysis was not limited to determining an average 
Po value. Instead, as in our previous studies |17|, re- 
sults for d, t, 3 He, and 4 He, were derived as a function of 
Wsurf, the velocity of the emerging particle at the nuclear 
surface, prior to Coulomb acceleration. From the rele- 
vant Pq values we then determined volumes using equa- 



FIG. 2: Equilibrium constants K c derived from the experi- 
mental yields. See text. Solid squares- K c (d), Solid triangles- 
K c (t), Open circles- K c ( 3 He) and Open squares K c (a). 
Solid diamonds indicate the values of K c (a) corresponding 
to intergral values of the temperature from 5 to 11 MeV (left 
to right) . The shaded bands indicate the uncertainity in den- 
sity. See text. 



tion ([2]). A comparison of these volumes indicated good 
agreement for t, 3 He and 4 He. The volumes derived from 
the deuteron data are typically somewhat smaller. This 
appears to reflect the fragility of the deuteron and its 
survival probability once formed 27| ■ For this reason we 
have used average volumes derived from the A > 2 clus- 
ters in the analysis. Given that mass is removed from 
the system during the evolution, we determined the rele- 
vant masses for each volume by assuming that the initial 
mass of the source was that obtained in the source fit- 
ting analysis and then deducing the mass remaining at 
a given v SU rf from the observed energy spectra. This is 
also an averaging process and ignores fluctuations. Once 
these masses were known they were used to determine 
an excluded volume for the particles. Addition of this 
excluded volume to the free volume produced the total 
volumes needed for the density calculations. These were 
determined by dividing the mass remaining by the total 
volume. This was done as a function of u SU rf • 
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RESULTS 
Temperatures and Densities 

Inspection of the results for the four different systems 
studied revealed that the temperatures, densities and 
equilibrium constants for all systems are the same within 
statistical uncertainties. Therefore we have combined 
them to determine the values reported in this paper. 

We present in Fig.[T]the experimentally derived density 
and temperature evolution of the IV source. Densities are 
expressed as total number of nucleons (including those 
in clusters) per fm 3 . There is a strong correlation of 
increasing temperature with increasing density. 

Estimated errors on these temperatures are 10% below 
p = 0.01 increasing to 15% at p ~ 0.03. The error in 
the derivation of the density arises from the uncertainty 
on the volume which is dominated by the uncertainty in 
temperature and the uncertainty in source mass derived 
from source fitting to complex spectra. The estimated 
errors are ±17%. We have included in each figure shaded 
bands representing this uncertainty in density. On the log 
plot the uncertainties in K are much smaller. 

Equilibrium Constants 

As already mentioned above, the absolute cluster 
yields and mass fractions calculated by the models de- 
pend upon the N/Z ratio of the matter, the model spe- 
cific nucleon-nucleon interaction assumed, and the vari- 
ous approximations of a given model. In addition, as all 
of the treatments assume chemical equilibrium, they also 
depend upon the number and type of competitive species 
included in the calculation. Historically, most EOS ta- 
bles, used in astrophysical simulations, have treated only 
neutronSjprotons, a-particles, and a single average heavy 
nucleus More recently some models have been de- 

veloped in which the number of competing species treated 
has been greatly expanded In an equilibrium situ- 

ation, all relevant equilibria must be simultaneously sat- 
isfied. Thus, if relevant species are not included, the 
calculated mass fractions of all included species will be 
in error. 

For this reason we do not believe that a direct compari- 
son with calculated mass fractions is the appropriate way 
to test the individual models. We choose rather to com- 
pare the experimentally derived equilibrium constants for 
the production of a particles with those of the models. 
The model derived equilibrium constants should be inde- 
pendent of proton fraction and the choice of competing 
species in a particular model. Specifically we define the 
equilibrium constant, K c , as 

KJA. (3) 
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FIG. 3: Comparison of experimental values of K c (a) with 
those from various EOS calculations. See text for the theo- 
retical models. The shaded band indicates the uncertainty in 
density for the experimental data. See text. 



where p(A,Z), p p and p n are respectively the densities 
of cluster of mass A and atomic number free protons 
and free neutrons. We express the densities in units of 
nucleons/fm 3 . In Fig. [5] we present the equilibrium con- 
stants, K c , for d, t, 3 He, and 4 He cluster formation as a 
function of density. For the purpose of later comparisons 
with EOS calculations we have interpolated the experi- 
mental data for the alpha cluster to determine the equi- 
librium constants at integral temperatures from 5 to 11 
MeV. These values are indicated by solid black diamonds 
on the figure. Temperatures increase as the density in- 
creases. Statistical errors are smaller than the symbols. 

For the purposes of a general assessment of the avail- 
able EOS calculations at low density, we present in Fig. [3] 
a comparison between the experimentally derived values 
of K c (a) and those predicted by available astrophysical 
EOS calculations. Note that this plot is logarithmic in 
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K c . The solid black diamonds connected by a dashed line 
indicate the experimental results. For a further compar- 
ison we show with a dotted line the results of the Virial 
EOS of Horowitz and Schwenk This EOS is based on 
virial coefficients derived from scattering data and has 
been suggested as a benchmark for other calculations at 
low density. Not surprisingly the calculated values of the 
equilibrium constant tend to converge at the lower den- 
sities. Even at the lowest densities sampled, however, 
there are significant differences. At higher densities, 0.01 
to 0.03 nuc/fm 3 , the various interactions employed all 
lead to a decrease of K c below that of the NSE values of 
Typel et al. as expected. However there are significant 
differences in the models, both in absolute value and in 
the rate of change of K c (a) at the higher densities. In 
the density range above ~ 0.01 nuc/fm 3 the predictions 
of the Virial model of Schwenk and Horowitz |5j, the 
Model of H. Shen et al. Q, the Lattimer Swesty Model 
with a Skyrme interaction having an incompressibility of 
375 MeV [3j and those of G. Shen et al. which employ 
the FSUGold [HI or NL3 interaction 9] all exhibit shal- 
lower slopes than the data and values all well above the 
data. 

The model predictions showing better agreement with 
the data above 0.01 nuc/fm 3 are the Lattimer and Swesty 
calculations employing Skyrme interactions with incom- 
pressibilities of 180 MeV and 220 MeV 3]. They are 
close to the data except for the two points for 180 MeV 
at highest density. 

The QSM model results are very close to the experi- 
mental results. The present calculations include the lat- 
est QSM estimates of momentum dependent in medium 
binding energy shifts @, Although well below the 
data and the other model calculations at the lowest veloc- 
ities, the results from reference Q, an NSE model which 
includes excluded volume effects are in good agreement 
with the data at the higher velocities. In a recent paper 
the extent to which the excluded volume concept can sim- 
ulate the in-medium effects, in particular Pauli blocking 
has now been explored [29] . 

Calculations using the recently improved RMF model, 
taking into account continuum effects 0, H^, produce 
predictions somewhat below the experimental data. In- 
corporation of the latest QSM model derived binding en- 
ergy shifts 28] into this model have not yet been made. 

While at the lowest densities the results of the the- 
oretical calculations converge towards the NSE result, 
the experimentally derived equilibrium constants at these 
lowest densities are significantly higher. This may in- 
dicate a limit in the experiment and inclusion of non- 
coalescence source of the alpha particles. Alternatively 
additional contributions to the alpha yield resulting from 
quartetting correlations which favor boson condensation, 
not treated in the models, might be indicated 31]. Fur- 



SUMMARY AND CONCLUSIONS 

Calculated equilibrium constants for specific cluster 
formation should be independent of n/p ratios and num- 
bers or types of other species included. Even so, reported 
model calculations of the equilibrium constant for alpha 
formation in low density nuclear matter vary by about 
two orders of magnitude in the density and tempera- 
ture region explored by the data presented in this paper. 
These new experimental data for these equilibrium con- 
stants provide important constraints on the low density 
equations of state. The use of the NSE which neglects 
in medium effects is clearly excluded as are several other 
models. While it is possible that these models might be 
improved , the data strongly indicate that accounting for 
in-medium effects, as in the semi-empirical excluded vol- 
ume approximation or the more sophisticated QSM 
approach [28l. l32jj . is required. This is particularly evi- 
denced in the density range of 0.01 to 0.03 nuc/fm 3 in 
the present data. 

We recall again that details of the low density compo- 
sition can be important in modeling supernova evolution 
and neutron star properties 0, 0, 33 1 . As emphasized in 
this paper, mistakes in K c (a) translate into errors not 
only in the yields of alpha particles but also in the yields 
of all competing species, regardless of the completeness 
with which those other species are included. The errors 
on the equilibrium yields of heavier species will increase 
with increasing mass. 
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ther investigations of this region might prove fruitful. 
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